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Power measurement is fundamental to transmission quality testing. 
This measurement need is extended to signals represented by digital 
bit streams. Accurate and precise measurements over a 60-dB range 
have been made using the digital signal processor. One algorithm 
that has been used measures the power of fixed- length sample se- 
quences. A second algorithm yields periodically updated power mea- 
surements of infinitely long sample sequences, but with slightly in- 
creased measurement ripple and frequency restriction. Theoretical 
expectations for measurement variation in the fixed- length measure- 
ments of noise power are also discussed. 

I. INTRODUCTION 

The measurement of power is fundamental to transmission quality 
testing. Power measurements of single tones, such as the milliwatt 
standard, are used to adjust transmission levels. Multiple-tone power 
measurements are used in nonlinear distortion testing. Examples of 
power measurements of band-limited noise are return-loss and C- 
message weighted noise measurements. 

This paper gives an overview, discusses the theoretical accuracy and 
precision of digital noise power measurements, and presents some 
results using the digital signal processor (dsp) A3990 for making power 
measurements. 

II. OVERVIEW 

The measurement of power will be presented following the block 
diagram of Fig. 1. Since the incoming signal is generally encoded for 
bit compression, the signal samples first must be decoded to linear 
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Fig. 1 — Time domain technique of power measurement. 

samples, as indicated in the first block of the figure. For ju-255 encoding 
and decoding, the dsp has a dedicated instruction set and associated 
circuitry. 

The decoded digital signal is scaled and passed through a digital 
filter. This paper discusses measurements with flat weighting, with C- 
message weighting, and through a C-notched filter. 

Note that Fig. 1 depicts essentially a time-domain approach, which 
should be contrasted with the frequency-domain approach depicted in 
Fig. 2. These two approaches are tied together by Parseval's Theorem: 



Signal Power = lim — 



y l {t)dt=\ S(w)dw = R(0), (1) 



where 



3~= averaging interval, 
y(t) = analog signal, 
S(w) = the power-density spectrum, and 

R(0) = autocorrelation function at zero lag, i.e., the dc component 
ofyHt). 

The first integral in eq. (1) is implemented as in Fig. 1; the second, as 
in Fig. 2. Implementation of the latter is not discussed in this paper. 

To compute the first integral, the sample values must be squared, as 
indicated in Fig. 1. Two different methods of determining the power 
from the squared signal were used. 

The first method is the straightforward approach. If the incoming 
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signal is considered stationary, the power can be approximated from 
the first integral of eq. (1) as follows: 



, r* -, n-i i n-i 

5(0)4 y 2 (t)dt^-±-ZylT = -Zyh 

y j Q iV i ,=o iv ,_ 



(2) 



where 

N = the number of samples, 

yi = ith signal sample, and 

T = interval between the samples (T= 125 /is in most voice 
telephone applications). 
This approach is termed block averaging (ba) in this paper. 

To use ba, N must be chosen large enough so that the measurement 
variation is within the required tolerance. In the next section, the 
probability of measuring noise power within certain confidence levels 
for different values of N is derived. 

The second method for extracting the dc component employs a 
convolution of the squared sample values with the infinite impulse 
response of a first-order, low-pass filter. This algorithm, described in 
Appendix B, has been termed exponentially mapped past (emp). 

The emp is not as applicable as ba because of the extra frequency 
components generated by squaring a signal. Because of the sampled 
nature of the signals, some components can be aliased into the pass- 
band of the emp low-pass filter and, thus, yield measurement ripples. 
For example, the emp parameters discussed in this paper yield a ripple 
of ±0.4 dB for a 15-Hz tone and a ripple of less than ±0.1 dB for an 80- 
Hz tone. 

Once the dc component has been extracted, it can be converted to 
a dB measurement before it is reported. A method for calculating the 
required logarithms with the dsp is described in Appendix C. 

Currently, the ba program reports a dB computation once per block, 
i.e., once every N samples, where N = 4096. However, when emp is 
used, linear-to-dB conversions can be made more frequently. After a 
conversion is made in the current emp program, the next conversion 
can be made after another three samples have been read from the 
input buffer. The rate with which conversions can be made and 
reported during emp power measurements depends on several vari- 
ables, as explained in Appendix D. 

For display or printout routines, a binary number representing a dB 
level can be converted to binary-coded decimal (bcd). A bcd routine 
was used to yield the ba signal-to-quantizing noise ratio measurements 
described below. This routine is not discussed in this paper. 

III. THEORETICAL PRECISION OF NOISE POWER MEASUREMENTS 

In this section the noise power measurement precision that is 
theoretically possible by digital power measurements is presented. 1 
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Consider the finite set of noise samples no, n u ... hn-i. The ac power 
of these N samples is expressed as 

i N-l 



where 



Pn(N) =j,Z (m - n)\ 



-i N-l 



(3) 



is the dc component of the noise. 

How large should N be in order that the estimated power P n (N) be 
within plus or minus some S n (in dB) of the actual power? That is, for 
a given N, what is the probability P{ } that 



,A1 P "M I" 1A1 PniM) 

10 log — hm 10 log — — 

ir ref M-.cc Xref 



<&> 



(4) 



where P ie { is any reference power? The following analysis assumes that 

the noise is Gaussian, but results are in most cases applicable to other 

types of noise (e.g., quantizing or coding noise). 

Assume the noise sample ra, to be an independent, zero-mean, 

Gaussian random variable with finite variance a 2 . Then P n (N) is an 

N — 1 
estimate of a 2 , with expected value = — — — a 2 (see Ref. 2). To the 

accuracy required for the ba program, P n (N) is effectively an unbiased 
estimate for N > 100. For a given 5„, eq. (4) becomes 



101og ^)- 101 o g ^ 

"ref "ref 



<8 n 



or 



CT 2 1() - 6n /10 < p n{N) < 2 1()5n /lO 



(5) 



(6) 



Now if the probability density function f n {a) of P n (N) can be found, 
then the probability P{ } of satisfying the inequality (5) is 



101og ^- 10 lo g ^ 
.r ref "ref 



•o 2 10*»/'° 

< 8„ } = I /"„(«) da. 

o 2 10"V 10 



(7) 



In Appendix A the density function is derived along with the proba- 
bility. The result is 



P{\~Vm-^a\<8 n } S 



h e - t 2/2 



dt, 



(8) 



where 
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Jfm = 10 log 



Pn(N) 



•> a = lOlog— , N>30, 

"ref 



h = y/2(10) s " /i0 N - V2JV-3, 



and 



/, = V2(10)- fi " /U) iV - V2N-3. 



For a graphical representation of eq. (8), see Fig. 3. For example, with 
500 samples the probability that the noise power measurement preci- 
sion is within ±0.2 dB is 53 percent. To meet standard specifications 
for noise and signal power measurements, 3 N was chosen to be 4096, 
yielding a precision of ±0.5 dB. 

IV. TEST RESULTS 

The ability to measure power precisely with both the ba and emp 
schemes can be seen by comparing ba and emp measurements with 
the actual signal-to-quantizing noise ratios (snrs) of ideally-encoded 
sine waves at levels ranging from 3 to —64 dBm. Such encoded sine 
waves have inherent quantizing noise frequency components across 
the voiceband range of frequencies. To make snr measurements, the 
C-notched and C-message filters are generally used. 
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Fig. 3 — Theoretical noise power measurement precision. 
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The C-notched filter is actually two cascaded filters, each imple- 
mented with three cascaded, second-order sections. The first filter, the 
C-message filter, has approximately unity gain from 1000 to 2500 Hz, 
and its attenuation increases gradually on either side of the passband 
to 54.7 dB at 60 Hz and to 13.7 dB at 3900 Hz. The second filter, the 
notch filter, is designed for attenuation of more than 50 dB from 989 
to 1020 Hz. The frequency response of the digital C-notched filter was 
measured by both the emp measurement routine and by an analog 
power meter. The results, limited by quantizing noise, appear in Fig. 
4. 

To determine the snr of an encoded signal between 1004 and 1020 
Hz, the C-message weighted measurement must be subtracted from 
the C-notched measurement. By means of a real-time developmental 
tool, the dspmate, snr measurements were made with the ba and emp 
programs using ideal, encoded sine waves at 1015.625 Hz. The ba and 
emp measurements yielded a range of ±0.5 dB, which was within the 
theoretical noise power measurement precision. In Fig. 5, the ba 
derived snrs are plotted against the actual snrs. 

In order to retain significant bits at low power levels, the programs 
were modified when the test signal powers were below -27 dBm. After 
/x-to-linear conversion, sample values from signals above —27 dBm in 
power were divided by 4, while sample values from signals below —27 
dBm in power were multiplied by 8. The dsp can be programmed to 
choose the appropriate scaling. 

Because of quantizing noise, the maximum snr is approximately 40 
dB. Thus, the maximum snr in Fig. 5 is comparable to the quantizing 
noise floor of Fig. 4. 
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V. CONCLUSION 

Two reliable methods for measuring power with the dsp have been 
presented: the ba approach, and the emp approach. Both of the 
approaches, when examined mathematically, have no significant bias 
in their expected values. The accuracy, therefore, is as good as the 
digital samples representing the signals being measured. 

The precision of signal power, noise power, and snr measurements 
were investigated. Digital signal processing measurements of these 
parameters showed that: 

(i) In general, signal power measurements were precise to within 
0.1 dB. 

(ii) With appropriate scaling, snr measurements were precise to 
within 0.5 dB over a 60-dB range. 

(Hi) Block-averaging noise power measurements all fell within the- 
oretical limits. 

In the emp program, frequencies in the range 80 to 3920 Hz yielded 
measurement ripples less than ±0.1 dB, and frequencies in the range 
30 to 3970 Hz yielded ripples less than ±0.2 dB. These measurements 
were made with an emp convolution having a 3-dB cutoff at 2.489 Hz. 
However, these ripples were about twice as large as the measurement 
ripples from the ba program. 
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Two advantages of emp over ba are compactness of code and ability 
to update a measurement 1365 times as often. Frequent updating may 
aid in identifying particular types of problems and, hence, aid in 
problem sectionalization. 



APPENDIX A 

Derivation of Noise Power Measurement Precision 

In the following analysis, the probability density function of P n {N) 
is presented and used to find, in computable form, the probability that 
the ac noise power P„(N) is within some 8„ (in dB) of the noise 
variance a 2 . 

Recall the definition of ac noise power: 

Pn(N) = i I (n, - n)\ 

N »=o 

As a result of assumptions that m is a zero-mean, Gaussian random 
variable with finite variance a 2 , the probability density function of 
P n (N) (see Ref. 4) is 

<JV-3)/2 -aNflo 2 

fn(a)= = , O) 

2 <N - 1,/2 (a/v^) A '- 1 r[(N-l)/2] 

which is a chi-square density. Hence, the probability P{ } that P n (N) 
is within 8„ of the actual noise power is 



101 o g ^>-10>o g ^ 

"ref " ref 



■o 2 10 6 " /l0 

< 8 n } - i f n (a) da. (10) 

'o*10-V 10 



Using the substitution t = aN/o 2 yields 

.2 



XOlogf^-iOlogf- 

"rc-f " ref 



<8 n 



ioV'on 



(N-3)/2 a -t/2 



dt 



io-« n /'°iv 



2 lN - 1)/2 T[(N - l)/2]* 



(11) 



For N > 30, this probability can be expressed in terms of the normal 
probability integral: 5 



"ref ^ref 



-t'/2 



<8 n \= | —;=du (12) 



/, v^ 



where 



h = V2(10) 5n/10 N - V2N-3 
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and 



lx = s/2(10)- s " /w N - V2AT - 3. 
Eq. (8) and Fig. 3 in the text follow. 



APPENDIX B 
Derivation of EMP 



The result of passing an analog signal y 2 (t) through an analog filter 
with impulse response g(t) is P(t), where 



Pit) = y\u)g(t - u) du. (13) 

Let T be the sampling period of the digital signal y 2 {nT), or y 2 (n) for 
short. The result of passing y 2 (n) through a digital filter with impulse 
response h{n), is P(n), where 

P(n)= I y 2 (k)h(n-k). (14) 

k CO 

An analog, first-order, low-pass filter uses 

... (Ae~ t/T , t>0 

git) = \o, t<0, < 15 > 

where t is the time constant of the filter. Such a filter has a 3-dB cutoff 
at (2t7t) _i Hz. 

If the impulse invariance h{n) = g(nT) is used to form an equivalent 
digital filter, the corresponding impulse response is 

,. , Ue~ nT/T , n>0 

h{n) '[0, n<0, (16) 

where T is the sampling period. Because of sampling, aliasing is 
introduced, but the effects of the aliased components are negligible if 
the cutoff frequency is low. 

For ease in notation, let m = T/r. From eqs. (14) and (16), 

P(n)= l y 2 (k)h(n-k), 

or 

n 

P(n) = £ y 2 (k)A e - m <»-*>. (17) 

From this, a simple recursive relationship for P(n) can be developed: 
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P(n + 1) - J /(*)A e — ^^ 

= e- m P(ra) + e _ V(» + D^ <? m (18) 



or 



P(n + 1) = e _m P(n) + Aj 2 (n + 1). (19) 

In eq. (19), A should be chosen to ensure unity gain at dc. To determine 
A, let 

' 2 <»> = {o,' "<o and «-»-«■ ,20> 

Then, 

P(n + 1) - e- m P(n) + AL, (21) 

which implies that P(n) is a geometric series: 



P(n) =AL%e 

k=0 



—mk 



,-m(n+l) 



-All^r-* <22) 

which approaches L for n approaching infinity if A = 1 - e~ m . 

The emp power measurement P(n) can, thus, be obtained by the 
recursion formula 

P( n ) = e - m P(n - 1) + (1 - e- n )y 2 (n), (23) 

where m = T/r. 

If the range of the signal power to be measured is large and no 
automatic gain control is to be used, then some double-precision 
arithmetic has to be done to save the least-significant bits resulting 
from sums and products. In particular, the result of squaring the input 
samples nearly doubles the number of significant bits in the accumu- 
lator. Therefore, all bits resulting from squaring should be saved. In 
addition, emp, which follows the squaring, should be implemented with 
double precision. 

This implementation is facilitated by representing e~"\ which is 
nearly unity for a low cutoff frequency, by 1 - 2" R , where R is an 
integer. The time constant t is 

t = -T/ln(l - 2~ R ). (24) 

For R = 9, t = 63.9 ms; and the 3-dB cutoff is 2.489 Hz, which yields 
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power measurements with ripples of less than ±0.1 dB at frequencies 
between 80 and 3920 Hz. 

To implement emp in double precision, the stored value P(n — 1) 
must occupy two storage locations. Conveniently, one location could 
contain its integral part, while the other, its fractional part. 

APPENDIX C 

Using the DSP to Compute the Natural Logarithm of a Number 

Using the dsp, the logarithm of up to a 40-bit number P can be 
computed to an accuracy of 0.001. Suppose that 

In P = ln[^(20] 

= \n.Jf+ S\n2, (25) 

where Jt is a real number and S is an integer. Then, In P can be 
computed from 

(i) a series expansion on Jt plus 
(ii) a table of multiples of ^ln 2. 
If < Jt < 2, the following expansion can be used: 

(Jt-1) 2 (Jt-lf (Jt-1) 4 
hiJt=(Jt-l)-± — -— - + - — — --- — — - + ••• (26) 

4U O T 

As indicated in Fig. 6, if 0.68 < Jt < 1.36, then In Jt can be computed 
to an accuracy of ±0.001 with only four terms. Since the upper bound 
on^is twice its lower bound, ^can be determined by repeated scalings 
of P by 2. 

APPENDIX D 

Rate of Linear-to-dB Conversions in the EMP Program 

In the current emp program, a linear-to-dB conversion of the power 
measurement is made after every third sample. Although every sample 
is used to update the linear power measurement, there is insufficient 
time to make a conversion to dB after every update. 

The number of samples, S, that must be used to update the linear 
power measurement S times and to make one linear-to-dB conversion 
can be determined from the following relation: 

S(R + U) + C < ST, (27) 

where 

R = the time to check whether the input buffer is full and, if so, 

to read it, 
U = the time to update the linear power measurement using the 

new sample, 
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C = the time to make a linear-to-dB conversion of the power 

measurement, and 
T = the sampling interval. 

Then S can be any integer greater than C/T - R - U.ln Fig. 7, the 
case for S = 3 is depicted. In this figure, L is the time necessary to load 
the input buffer. 

As shown in Fig. 7, the second sample in each group of three must 
be read before the third sample begins to be loaded. Expressed sym- 
bolically, this means 
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R+ U+C + R < T+T- L, 



or 



C < 2T - L - 2R - U. 



(28) 



If eq. (28) were not satisfied by C, then the second sample in each 
group of three would have to be stored during the linear-to-dB con- 
version and then the sample value returned after the conversion was 
complete. 

In the current emp program, the upper bound on C was sufficiently 
high that, with S = 3, each linear-to-dB conversion could be followed 
by a conversion to BCD. 
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